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Statistical  signal  processing  functions  such  as  signal  detection,  estimation,  and  identification  play 
a  key  role  in  the  development  of  effective  communications,  radar,  and  sonar  systems.  For  example, 
advanced  statistical  methods  are  emerging  as  being  particularly  important  in  digital  communications 
systems  operating  in  channels  corrupted  by  interference  from  such  phenomena  as  multiple-access  noise, 
intentional  jamming,  and  impulsive  noise  sources.  Conventional  demodulation  methods,  such  as 
coherent  matched  filtering,  often  suffer  serious  performance  degradation  when  subjected  to  interference 
of  these  types;  however,  this  degradation  can  frequently  be  eliminated  through  the  use  of  more  sophis¬ 
ticated  signal  processing  techniques. 

A  central  issue  in  the  design  of  effective  signal  processing  procedures  for  systems  operating  in 
channels  such  as  those  noted  above  is  that  of  channel  identification.  Although  certain  aspects  of  chan¬ 
nel  identification  have  been  studied  extensively,  one  area  in  which  there  is  a  pressing  need  for  further 
research  is  that  of  identification  of  impulsive  channels.  Communication  systems  are  seldom  interfered 
with  by  white  Gaussian  noise  alone,  yet  receiving  systems  in  general  use  are  those  which  are  optimum 
for  white  Gaussian  noise.  The  man-made  electromagnetic  environment,  and  much  of  the  natural  one, 
is  basically  "impulsive,"  i.c.,  it  has  a  highly  structured  form  characterized  by  significant  probabilities  of 
large  interference  levels.  In  addition  to  the  man-made  electromagnetic  environment,  there  are  many 
other  different,  common  and  widely-used  communications  and  remote -sensing  type  channels  where 
impulsive  noise  dominates,  c.g.,  cxtra-low-frcqucncy  (ELF)  channels,  urban  radio  networks,  underwater 
acoustic  channels,  and  telephone  line  channels.  This  is  in  contrast  to  the  Gaussian  noise  processes 
inherent  in  traasmitting  and  receiving  elements.  Since  the  conventional  receivers  arc  effectively  linear, 
the  impulsive  character  of  the  interference  can  drastically  the  performance  of  conventional  systems. 
Furthermore,  since  it  has  been  well-established  that  the  performance  of  communications,  radar,  and 
sonar  systems  can  be  greatly  enhanced  if  the  true  statistics  of  the  channel  arc  known  and  exploited,  the 
problem  of  identification  of  impulsive  noise  channels  is  important  in  the  development  of  systems  that 
can  approach  the  performance  limits  set  by  such  channels. 

This  research  project  is  concerned  with  the  design  and  development  of  new  nonlinear  algorithms 
for  identification  of  nongaussian  noise  processes  and  for  reception  of  signals  corrupted  by  nongaussian 
noise.  In  particular,  the  overall  objectives  of  this  research  project  arc;  (i)  the  development  of  suitable 
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channel  models  for  impulsive  interference;  (ii)  the  derivation  and  analysis  of  new  globally  optimum 
nonlinear  estimation  procedures  for  the  identification  of  impulsive  interference  channels,  and  (iii)  appli¬ 
cation  of  the  developed  algorithms  in  the  design  of  optimum  and  sub-optimum  adaptive  receivers  for 
the  detection  of  signals  in  non-Gaussian  noise  environments. 


RESEARCH  SUMMARY 


During  this  reporting  period,  the  focus  of  our  work  has  been  on  the  problem  of  obtaining 
optimum  and  efficient  identification  and  detection  procedures  for  impulsive  channels.  Of  particular 
interest  is  the  Middleton  Class  A  noise  model,  which  is  a  widely-accepted  statistical-physical  model  for 
impulsive  interference  superimposed  on  a  Gaussian  background.  The  model  has  two  basic  parameters 
that  can  be  adjusted  to  fit  a  wide  variety  of  impulsive  noise  phenomena  occurring  in  practice. 

We  have  considered  the  problem  of  basic  batch  estimation  of  the  Class  A  parameters.  In  particu¬ 
lar.  we  have  developed  a  consistent  and  asymptotically  efficient  estimator  of  the  Gass  A  parameters 
that  is  practically  efficient  as  well  [17], [18].  The  development  of  this  estimator  began  with  an 
improved  method-of-moments  estimator  of  the  parameters  based  on  the  extraction  of  two  moments 
chosen  for  their  sensitivity  to  the  channel  parameters.  We  have  shown  that  this  moments  estimator  is 
relatively  simple  to  implement  and  yields  unique  estimates  of  the  parameters  that  converge  to  the  true 
parameters  at  near-optimal  rates.  Wc  have  also  shown,  via  an  extensive  simulation  study,  that  the  esti¬ 
mator  yields  good  estimates  of  the  parameters  for  moderate  sample  sizes  but  that  a  significant  improve¬ 
ment  in  performance  over  the  method-of-moments  estimator  is  possible  if  an  asymptotically  efficient 
estimator  can  be  found.  The  search  for  such  an  estimator  began  with  a  search  for  a  method  which 
identifies  the  consistent,  and  hence,  asymptotically  efficient  sequence  of  roots  to  the  likelihood  equa¬ 
tion.  One  method  that  has  the  potential  of  providing  such  an  estimator  is  maximum  likelihood.  Wc 
have  shown  in  this  study  that  the  maximum  likelihood  estimate  (MLE)  is  consistent  for  the  Gass  esti¬ 
mation  problem,  despite  the  existence  of  multiple  roots  in  the  asymptotic  likelihood  equation  [6j. 
However,  the  existence  of  these  multiple  roots  makes  the  MLE  very  difficult  to  implement  For¬ 
tunately,  the  consistent  sequence  of  roots  can  be  located  using  an  alternative  method.  This  method 
consists  of  initiating  Newton’s  root-finding  method  on  the  likelihood  equation  wfith  the  method -of- 
moments  estimator  and,  in  so  doing,  both  the  consistency  of  the  moments  estimator  and  the  efficiency 
of  likelihood-based  estimation  arc  retained.  We  have  shown  that  this  Momcnt/Likclihood  procedure 
does  indeed  yield  consistent  and  asymptotically  efficient  estimates  of  the  parameters.  Moreover,  via  an 
extensive  simulation  study,  it  was  seen  that  the  proposed  Moment/Likelihood  procedure  yields  excel¬ 
lent  estimates  of  the  parameters  for  moderate  sample  sizes.  In  addition,  wc  have  developed  an  optimi¬ 
zation  scheme  based  on  the  simplex  method  that  can  be  shown  to  converge  to  a  relative  maximum  of 
the  likelihood  function.  An  extensive  simulation  study  of  this  scheme  indicates  that  it  performs  very 
well  for  small  sample  sizes. 

The  problem  of  recursive  estimation  of  the  Class  A  parameters  was  also  investigated  [7],  Our 
primary  objective  in  this  phase  of  the  study  was  the  development  of  recursive  channel  estimators  that 
arc  efficient  for  small  sample  sizes.  Only  in  this  way  can  wc  guarantee  that  our  channel  estimator  is 
capable  of  tracking  rapidly-varying  nongaussian  channels.  In  particular,  wc  have  developed  a  global 
recursive  estimator  of  the  parameters  that  can  overcome  the  shortcomings  of  conventional  recursive 
stochastic  gradient-type  algorithms  which  make  them  unsuitable  for  impulsive  channel  identification. 
The  primary  shortcoming  is  that  the  conventional  algorithms  cannot  uniquely  identify  the  Gass  A 
parameters  due  to  the  existence  of  multiple  roots  in  the  asymptotic  likelihood  equation.  This  recursive 
estimator  is  based  on  the  use  of  a  nonstandard  stochastic  approximation  procedure  wherein  the  standard 
Fisher’s  information  matrix,  which  governs  the  step  size  in  the  recursion,  is  replaced  with  a 
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"complete-data"  information  matrix.  Two  initiating  estimators  for  this  complete-data  stochastic  approx¬ 
imation  procedure  were  developed.  The  first  is  a  hatch  estimator  which  extracts  estimates  of  the 
parameters  directly  from  a  histogram  of  the  Class  A  envelope  data.  The  second  is  a  melhod-of- 
moments  estimator  applied  to  a  nontraditional  set  of  observations,  referred  to  as  "complete-data."  At 
each  iteration  of  this  initiating  recursive  estimator,  complete  data  are  generated  from  the  traditional,  or 
incomplete,  observations  using  the  estimate  of  the  parameters  obtained  at  the  previous  iteration.  A 
fixed  parameter  vector  in  the  parameter  space  initiates  the  process.  By  equating  sample  moments  of 
the  complete  data  to  the  true  moments,  estimates  of  the  parameters  are  then  obtained.  Both  initiating 
estimators  were  found  to  yield  estimates  which  lie  in  the  region  of  convergence  of  the  stochastic 
approximation  procedure.  Moreover,  via  an  extensive  simulation  study,  it  was  seen  that  the  proposed 
stochastic-approximation-based  estimators  yielded  global  recursive  procedures  through  which  excellent 
estimates  of  the  parameters  were  obtained  for  small  sample  sizes. 

During  this  period,  we  have  also  developed  and  analyzed  the  performance  of  several  rton- 
paramctric  density  estimators  and  detectors  for  impulsive  noise  channels  14), [5],  These  nonparametric 
density  estimators  exploit  the  known  properties  of  the  given  noise  environments  and.  in  so  doing,  yield 
estimates  that  closely  approximate  the  true  densities  for  small  sample  sizes.  Two  estimators  have  been 
proposed:  the  first  is  particular  to  envelope  densities,  while  the  second  is  designed  for  symmetric- 
instantaneous  amplitude  densities.  The  standard  automatic  kernel  and  variable  kernel  estimators  pro¬ 
vided  the  basis  for  the  construction  of  these  two  estimators.  By  incorporating  the  two-component 
structure  of  the  impulsive  noise  densities  to  be  estimated  into  the  framework  of  these  estimators  and 
choosing  the  variables  of  the  estimators  appropriately,  tremendous  improvements  in  performance  over 
the  basic  kernel  estimators  were  achieved.  In  addition,  a  method  for  evaluating  the  L  r performance  of 
an  "optimum"  parametric-based  density  estimator  was  derived.  Upon  using  this  method,  it  was  seen 
that  the  performance  of  the  proposed  nonparametric  density  estimators  was  comparable  to  that  of  their 
optimum  counterparts  for  small  values  of  the  sample  size  for  the  wide  variety  of  impulsive  noise  densi¬ 
ties  considered  in  this  study.  These  nonparametric  density  estimators  were  then  used  in  the  develop¬ 
ment  of  various  detection  strategies.  In  particular,  two  likelihood  ratio  tests  were  formulated  using 
these  density  estimators,  one  based  on  the  classical  robust  test,  the  other  on  the  /,  j -error-based  detector. 
Via  an  extensive  simulation  study,  it  was  seen  that  the  performance  of  these  schemes  was  very  near 
that  of  the  optimum  likelihood  ratio  test  for  small  values  of  the  sample  size.  The  proposed  non¬ 
parametric  density  estimators  were  then  used  in  the  development  of  a  weak-signal  detector  that  approx¬ 
imates  the  locally  optimum  nonlinearity.  By  using  the  basic  features  of  the  robust  and  L\ -error-based 
tests,  a  "prcccnsoring"  method  was  developed  for  this  detector  which  suppresses  the  effects  of  errors  in 
the  density  estimates.  Via  an  extensive  small-sample  size  performance  analysis,  the  performance  of  the 
resulting  detector  was  found  to  be  very  near  that  of  the  optimum  nonlinearity  for  the  wide  class  of 
densities  considered  in  this  study.  Furthermore,  for  a  variety  of  impulsive  noise  densities  this  detection 
scheme  was  shown  to  outperform  several  classical  nonparametric  detectors. 

A  nonparametric  density  estimator  for  positive  random  variables  was  also  developed  and  analyzed 
[8],  This  kernel-based  estimator  makes  use  of  a  basic  feature  of  positive  random  variables  in  its 
design,  namely  that  a  positive  random  variable  can  be  represented  as  the  norm  of  a  random  vector  hav¬ 
ing  arbitrary  direction.  Among  other  applications,  this  approach  is  well-suited  for  estimating  densities 
of  the  envelope  of  electronic  and  acoustic  waveforms  (c.g.,  the  norm  of  the  in-phase  and  quadrature 
components  of  signals).  Depending  on  the  dimension  of  the  normed  vector  space  dial  is  employed, 
this  estimator  also  assumes  different  forms,  each  form  being  applicable  to  a  different  class  of  densities. 
The  proposed  estimator  was  analyzed  extensively.  First,  the  estimator  was  shown  to  possess  some  of 
the  usual  desirable  properties  of  nonparametric  estimators,  such  as  pointwisc  and  L,-consistcncy. 
Then,  an  upper  bound  on  the  expected  value  of  the  L  j  error  was  derived  and,  using  a  minimax  stra¬ 
tegy,  the  optimal  kernel  and  smoothing  factor  for  the  estimator  were  determined.  In  addition,  the 
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appropriate  forms  of  the  estimator  for  the  different  classes  of  deasities  were  specified.  Finally,  the 
asymptotic  performance  of  the  proposed  estimator  and  the  standard  kernel  estimator  (restricted  to  die 
nonnegative  real  axis)  were  compared  using  the  derived  upper  bound  on  the  expected  value  of  the  /, , 
error.  Via  this  comparison,  the  proposed  scheme  was  seen  to  outperform  die  standard  estimator  for  a 
wide  variety  of  densities. 
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